Acta Phys. Hung. A 19/1 (2004) 000-000 

HEAVY ION 
PHYSICS 



in 
o 
o 

(N 

c 
3 



Ellipsoidal flows in relativistic hydrodynamics of finite 
systems 



Yu.M. Sinyukov 1 and Iu.A.Karpenko 1 ' 2 

1 Bogolyubov Institute for Theoretical Physics, Kiev 03143, Metrologichna 14b, 
Ukraine 

2 National Taras Shevchenko University of Kyiv, Kiev 01033, Volodymyrska 64, 
Ukraine. 



' Received 1 January 2004 

(N 

Abstract. A new class of 3D anisotropic analytic solutions of relativistic hydro- 
dynamics with constant pressure is found. We analyse, in particular, solutions 
' corresponding to ellipsoidally symmetric expansion of finite systems into vac- 

, mim. They can be utilized for relativistic description of the system evolution 

' in thermodynamic region near the softest point and at the final stage of the 

hydrodynamic expansion in A+A collisions. The solutions can be used also 
for testing of numerical hydrodynamic codes solving relativistic hydrodynamic 
equations for anisotropic expansion of finite systems. 
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H ' 1. Introduction 

The solutions of relativistic hydrodynamics are used typically in cosmology, astro- 
physics and for an analysis of the processes of ultra-relativistic heavy ion collisions. 
It is remarkable fact that in these so different fields the Hubble-like solutions play 
the important role. As for A+A collisions the Hubble-like models becomes to be 
popular for a description of the experimental data only lately [^| since it was found 
in RHIC experiments that rapidity distribution is not flat even at these top en- 
ergies, as it was expected (see, e.g., rewiev in [H]). Thus, one can suppose that 
the initial state could be more close to spherically symmetric one rather then to 
boost-invariant in longitudinal direction [EH E]- 

The spherically symmetric hydrodynamic solutions with the Hubble velocity 
distribution, v = r/t, has been considered for the first time in Ref. The equation 
of state (EoS) was chosen as ultrarelativistic one: p = c\e. Some generalization of 
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these results was proposed in a case of the Hubble flow for EoS of massive gas with 
conserved particle number in Ref. [ 6 . 

It is naturally, however, that, unlike to the Hubble type flows, the velocity 
gradients in thermal matter formed in A+A collisions should be different in differ- 
ent directions since there is an initial asymmetry between longitudinal and trans- 
verse directions in central collisions and, in addition, between in-plane and off-plane 
transverse ones in non-central collisions. The another important remark is that the 
Hubble-like hydrodynamic solutions is related to infinite systems while the matter 
in A+A collisions is occupied essentially finite region. 

In this letter we search for analytical solutions of the hydrodynamic equations 
describing 3D asymmetric relativistic expansions of finite systems. 

2. General analysis 

Let us start from the equations of relativistic hydrodynamics: 

d v T^ = 0, (1) 

where the energy-momentum tensor corresponds to a perfect fluid: 

T^ = (s+p)u fi u u -p-g fiv (2) 

We can attempt to find a particular class of solutions and therefore have to make 
some simplifications of 

Let us put pressure to be constant, 

V = Po = const, (3) 

in other words, we are searching for the solution near the softest point where the 
velocity of sound c s — 0. 

Then we find from Eq.QJ 

( f i P )a„(//) + iiVi ( = o (4) 

Contracting the latter equation with we obtain: 

(s + p)d v u v + u v d v e = 0. (5) 

and, thus, the remaining equation to satisfy is 

xFdvU* = 0, (6) 

which means that flow is accelerationless in the rest systems of each fluid element; 
this property holds for the known Bjorken (boost-invariant) and Hubble flows. 

Finally, the relativistic hydrodynamics at the softest point is described by the 
equations (|SJ and © for the hydrodynamic velocities u^, and energy density. A 
serious problem is, however, to find non-trivial solutions for the field ^(x) of hy- 
drodynamic 4- velocities. 
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3. Relativistic ellipsoidal solutions 

First the solution at p = const was obtained in [[7] as physically corresponding to 
a thermodynamic state of the system in the softest point with the velocity of sound 
c 2 = 0. Such a state could be associated with the first order phase transition. In 
A+A collisions it corresponds, probably, to transition between hadron and quark- 
gluon matter at SPS energies. The solution proposed in [ Ej has the cylindrical 
symmetry in the transverse plane and the longitudinal boost invariance: 

u fl = l(-,v-, v-,-), (7) 
t r r t 



where r = \Jt 2 — z 2 , 7 = (1 — v 2 ) 1 ^ 2 , r is transverse radius, r — y 1 x 2 + y 2 , and 



(8) 



1 + OLT 

describes axially symmetric transverse flow. 

The above solution has, however, a limited region of applicability since the 
boost invariance is not expected at SPS energies and can be used only in a small 
mid-rapidity interval [0, it is not reached even at RHIC energies [0- But most 
important is that in non-central collisions there is no axial symmetry and, therefore, 
one needs in transversely asymmetric solutions to describe the elliptic flows in these 
collisions, e.g., 1*2 coefficients. Now we propose a new class of analytic solutions of 
the relativistic hydrodynamics for 3D asymmetric flows. 

Firstly we construct the ansatz for normalized 4-velocity: 

„ = r * _ °fc(*)sfc _ X (Q) 

V* 2 -Eaif(*KV* a -Eaif(*K ? 

where the Latin indexes denote spatial coordinates, which are functions of time 
only. In this case a set of nonequal Oj induces 3D elliptic flow with velocities 
vi = di(t)xi/t: at any time t the absolute value of velocity is constant, v 2 = const, 
at the ellipsoidal surface ^l x \ = const. 

The condition © of accelerationless thus reduces to the ordinary differential 
equation (ODE) for the functions tii(t): 

dm _ gj - af 

dt t ' ( ' 



the general solution of which is: 



a^t) = * (11) 



where Tj is some set of 3 parameters (integration constants) having the dimension 
of time. The different values T\ , T2 and T3 results in anisotropic 3D expansion with 
the ellipsoidal flows. 
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The solution of Eq. (JSJ for energy density e is found as follows. Taking into 
account that d^u^ — y^Qj/T, where 



t 2 (12) 



one can get 

(e + ft,) ^ Oi(i) + ifte + ^ a t {t)x l d t e = 0. (13) 
General solution of the equation is 

where F e is an arbitrary function of its variables. At fixed Tj that define the 
velocity profile, function i 7 ^ is completely determined by the initial conditions for 
the enthalpy density. 

If some value, associated with a quantum number or with particle number in 
a case of chemically frozen evolution is conserved [ |S] , then one should add the 
corresponding equation to the basic ones. Such an equation has the standard form 

nd v u v + u^d^n = Q (15) 

where n is associated with density of the correspondent conserved value, e.g., with 
the baryon or particle density. A general structure of this equation is similar to Eq. 
((T3^l and, therefore, the solution: 



P ( x l x ? Z3 \ 

J "Vt+Ti' t+T 2 ' t+Tz> 

(t + TiXi + TaXt + Ta)' 



(16) 



which looks like as i|14|) and the function F n is also an arbitrary function of its 
arguments and can be fixed by the initial conditions for (particle) density n. The 
corresponding chemical potential is not zero and describes the deviation from chem- 
ical equilibrium in relativistic systems. The thermodynamic identities lead to the 
following expression for the temperature: 

where Ft is some function of its arguments that is defined by the initial conditions 
for e and n as well as by EoS e — e(n,T). If the initial enthalpy density profile 
is proportional to the particle density profile, F n (^, f|, f|) ~ F e (^, ||), then 
T = const and so n = const. 

In a case of chemically equilibrium expansion of the ultrarelativistic gas when 
the particle number is uncertain and is defined by the conditions and parameters of 
the thermodynamic equilibrium, e.g., by the temperature T, the chemical potential 
fi = 0. Then it follows that T — const for such a system and the entropy s = 
(e(t, x) +p )/T where e(t,x) is defined by ifHfl. 
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4. Relativistic anisotropic expansion of finite systems 

Let us describe some important particular solutions of the equations for relativistic 
ellipsoidal flows. If one defines the initial conditions on the hypersurface of constant 
time, say t = 0, then t is a natural parameter of the evolution. Such a representation 
of the solutions, which are similar to the Bjorken and Hubble ones, with velocity 
field Vi = aiXi/t has property of an infinite velocity increase at x — > oo. A real 
fluid, therefore, can occupy only the space-time region where |v| < 1, or f 2 > 0. To 
guarantee the energy-momentum conservation of the system during the evolution, 
all thermodynamic densities have to be zero at the boundary of the physical region, 
otherwise one should consider the boundary as the massive shell [ . Hence in 
the standard hydrodynamic approach the enthalpy and particle density must be 
zero at the surface defined by |v(i, x)| = 1 at any time t. One of a simple form of 
such a solution (for the case of particle number conservation) can be obtained from 

(11311,(0 choosing F e ^ n — cxp ( -fe 2 |^ ] : 



where f is defined by Ijl2|l . and the constants C £ , C n , b e and b n are determined 
by the initial conditions as described in the previous section. As one can see, 
the enthalpy density tends to zero when |x| becomes fairly large approaching the 
boundary surface defined by |v(i, x)| = 1, in the other words, when f — > 0. Thus 
the physically inconsistent situation when massive fluid elements move with the 
velocity of light at the surface f = is avoided. Of course, in such a solution one 
has to put a constant pressure to be zero, po = 0. 

On the other hand, the solution (|18|) , (|19|) can describe effectively finite systems 
at the hypersurface f = const. To show this one can substitute by the definition 
t 2 = f 2 + a i x i m the solution, and get: 



supposing that x% are spatial coordinates on the hypersurface. This describes effec- 
tively finite hydrodynamic flow with spatial radii Ri r/b e at t 3> TV 

As it follows from an analysis of the behavior of the thermodynamic values 
in the previous section, the temperature is constant if b e = b n , otherwise one can 
choose the temperature approaching zero at the system's boundary, e.g., for EoS 




(19) 



(18) 




(20) 



(21) 
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which is linear in temperature, the latter has the form 

T = const, b e = b n , 

(22) 

T ~ e ~ (b '~ fc " ) ^ -» 0, \v{x)\ -» 1, 6 e > &«• 



Note that in the region of non-relativistic velocities, w 2 = X) ""t^ 1 ^ 1) the 
space distributions of the thermodynamical quantities l(TB|) . ifTHj) has the Gaussian 
profile: 



^e"^^, (23) 



The forms of the solutions (|23|l are similar to what was found in Ref. [ 1TJ as the 
elliptic solutions of the non-relativistic hydrodynamics equations. In this sense the 
solution proposed could be considered as the generalization (at vanishing pressure) 
of the corresponding non-relativistic solutions allowing one to describe relativistic 
expansion of the finite system into vacuum. 

One can note that the case of equal flow parameters Tj = and b e = b n = 
induces formally Hubble-like velocity profile. 

It is worthily to emphasize that the physical solutions with non-zero constant 
pressure have a limited region of applicability at least in time-like direction: if one 
wants to continue the solutions to asymptotically large times, then (e +Po)t-^oo ~ 
S — ► 0, and this results in non-physical asymptotical behavior e — ► —po, unless we 
set po = 0. Therefore, it is naturally to utilize such kind of solutions in a region of 
the first order phase transition, characterized by the constant temperature and soft 
EoS, c 2 s — dp/de w 0, or at the final stage of the evolution that always corresponds 
to the quasi-inertial flows. 



5. Conclusions 

A general analysis of the flows at soft EoS p = const in the relativistic hydrody- 
namics is done. A new class of analytic solutions for 3D relativistic expansion with 
anisotropic flows is found. These solutions can describe the relativistic expansion of 
the finite systems into vacuum. They can be utilized for a description of the matter 
evolution in central and non-central ultra- relativistic heavy ion collisions, especially 
during deconfinement phase transition and the final stage of evolution of hadron 
systems. Also, the solutions can serve as a test for numerical codes describing 3D 
asymmetric flows in the relativistic hydrodynamics. 
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